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Abstract 

Cosmic Microwave Background (CMB) Anisotropics is a subject of intensive research 
in several fields of sciences jHD02j . In this paper we start a systematic development of 
basic notions and theory in statistics according to the application for CMB. The main 
result of this paper is the necessary and sufficient condition for isotropy of a non-Gaussian 
field in terms of spectra. Clear formulae between bi-, tri- and polyspectra and bi-, tri-, 
and higher order covariances are also given. 

Keywords: Bispectrum, Trispectrum, Angular poly-Spectra, Cosmic microwave back- 
ground radiation; Gaussianity; spherical random fields 

1 Introduction 

" Cosmological Principle (first coined by Einstein): the Universe is, in the large, homogeneous 
and isotropic" |Bar99j 

In the last decade or so there has been some growing interest in studying the spatial - time 
data measured on the surface of a sphere. These data includes cosmic microwave background 
(CMB) anisotropics [KKOGj . |()Hn2j . |AH12j . medical imaging |KBID93) |Kakl2j . global and 
land-based temperature data | Jon94j . |SRT06j . gravitational and geomagnetic data, climate 
model |N(]81j . 

One of the problem in focus is the Non-Gaussianity of the observations, which leads 
to the investigation of higher order angular spectra called polyspectra, [HuOlj . [BSMRIO] . 
|BLSH12] . Angular polyspectra, in particular the bispectrum and trispectrum shows to be 
an appropriate measure of Gaussianity since for a Gaussian process all higher (then second 
order) spectra are zero, |KSH11| . An other important question is the Monte Carlo simulation 
of non-Gaussian isotropic maps with a given power spectrum and bispectrum and possible 
polyspectrum, [(MOl] . [R,HS+n5] . 



In this paper some general properties of the angular polyspectra will be given for isotropic 
stochastic processes on sphere. Our treatment follows the basic theory of higher order spectra 
for non-Gaussian time series, |Bri65j . [BriOl] . [SRG84] . [Ter99j . The symmetry relations and 
the appropriate series expansion of cumulants might have influence on the estimation of 
polyspectra as well, |SK12] . |IMK12] . 
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The isotropy assumption implies very particular form for the angular spectra therefore a 
delicate question is the construction of isotropic stochastic maps on the sphere. 

2 Gaussian isotropic fields 

Gaussian isotropic processes on sphere has a long history starting with |Obu47] . some basic 
theory and references can be found in |Yag61| , | Jon63j . |McL86j |Yad83] |Yag87| . Due to 
the expansive recent applications there are new books [GGlOj . [CWll] and several papers 
covering a number of problems in general for spatial-time observations. 

We consider a stochastic process X (L) on the unit sphere §2 in I^^i where L = (??, 93), the 
co-latitude -i? G [0,7r] and the longitude (p G [0, 27r]. Let us suppose that X {L) is continuous 
(in mean square sense), then it has a series expansion in terms of spherical harmonics 

00 I 

1=0 m=-l 

where the coefficients are given by 

= f X (L) {L)n (dL) , 

where J7 (dL) = sin ■dd'ddip is Lebesque element of surface area on S2. It is generally ac- 
cepted, in time series analysis that X (L) called linear if ZJP are independent and identically 
distributed {EX (L) = 0). Now it is straightforward that from the linearity here follows 
that the covariance function C (Li, L2) = EX (Li) X (L2) depends on the angular distance 7 
between Li and L2 only. Necessarily the covariance function depends on the central angle 
between locations and have the form 

C (Li , L2) = C (cos 7) = 5; ^'^r^' ' 

£=0 

where Pi denotes the Legendre polynomial, see |Yad83] . and all coefficients fe > 0. Note 

here that the inner product Li • L2 = cos 7. In other words C (Li,L2) is invariant under 

the group of rotations, i.e. X (L) is isotropic. Moreover the assumption of linearity, i.e. the 

independence of the triangular array so strong that the Gaussianity also follows. In other 

words the only linear field on the sphere is the Gaussian one. We shall consider a linearity 

which looks weaker although it will be shown to be equivalent to the whole independence of 
'7m 

■ 

Definition 1 The field X [L] is linear if the generating array is uncorrelated and for 
fixed degree £, { Zp \ m = —£, —£ + 1 . . . ,i — 1,1} are independent. 

This concept of linearity corresponds to the physical approach to the theory of angu- 
lar momentum. One may consider linearity in some really weaker sense, namely the rows 
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of the array Z^, are independent but the variables Z"^ inside the row I are 
uncorralated. 

The convergence of the series^^Q /^^f^ is equivalent to the continuity of C (•) on [—1, 1]. 
The superposition ^ corresponds to the superposition of a covariance function on the real 
line in terms of its spectrum and the orthogonal system {exp {i2TrXk) , A; = 0, 1,2...}, hence 
we treat i as the frequency and Ci = fi is the value of the spectrum at £. Since Pi (cos 0) = 1, 
the variance EX (L)^ is broken up into the sum of spectra therefore we have the analysis of 
variance. The inversion 

Cf = C (cos 7) Pe (cos 7) Q (dL) , 

is also valid. The stochastic process X (L) itself has spectral representation ([1]) also in terms 
of spherical harmonics with complex values, detailed account for spherical harmonics Y^™ 
is found in |VMK88] . |SW71] . For a given X (L) we have the inversion 

Z^= [ X (L) Yp {L)n (dL) . 

The orthogonal random 'measure' is a triangular array for each &x i; m = —£, —i + 
1 . . . — i.e. rows contain 2i + l, i.i.d Gaussian random variables, EZ^^ = 0, EZ^^Z^* = 
fiS£,kSm,n- The reason that fi does not depend on m, is the Funk-Hecke formula [23l see 
|Yad83j .' |.Ton63j . 

The C (Li • L2) is strictly positive definite if all fi is > 0, and only finitely many are zero 
f [Sdl97] . [SM] ). 

Example 2 Matern Class of Covariance Functions \Mat86^ . ISGOSf Suppose one has a co- 
variance function on the real line then its restriction to [0, vr] gives a covariance function 
Cq (7) , on the sphere §2; Cq (7) = C (cos 7), for instance 

Co (7) = ( y) ^"^^^ , ^ > 0, 1. > 0, 

fi = J Co (7) (^°^ ^) '^'y- 

Where is the modified Bessel function of the second kind. Here the smoothness parameter 
V controls the continuity, as well as i) controls the regularity IGGIO^ . IShell]/ . Note Ky (r) ~ 
r {v) {r /2)~^ /2 if r ^ 0, and u > 0, hence Cq (7) /o"^ is a correlation function. 

Example 3 The generating function of Legendre polynomial Pi is 

00 

J2Piix)z^ = (1 -2xz + z^)"^^^ X G (-1,1) , \z\<l 
£=0 

see EH < 2: < 1, hence 

C (cos 7) = =, 
y 1 — 2zcos7 + 
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is covariance function on [—1,1], with fi = 21^^ ■ Similarly for n> 2 we have a covariance 
function 

C (cos 7) = o.(n-2)/2 ' 

(1-2ZC0S7 + Z2)(" 2)/2 

for any < 2; < 1. 

Example 4 From the well known series 

see JYadSS^ . we find the covariance function on [—1, 1] 

C (cos 7) = 

(l-2acos7 + a2)3/2 

for any < a < 1, and fi = Aira^. Similarly for n >2, we have 

C (cos 7) = —p^. 

(1 — 2a cos 7 + a^j ' 



Moreover the series 

°° t^ 2 

EPe (x) — = In =, \t\ < 1, \x\ < 1 

^ ^ i 1 - tx + Vi^^2^^tT^ ' ' ' ' 

see lPBM86\j . 5.10.1.4- Tom 2, gives a covariance function as well. 

Example 5 It is known, \Bre6S^ . ] WCOlf . that the probability density on the sphere 

C (cos 7) = , , , exp (k cos 7) , k > 

47rsmh(K) 

has series expansion 

where \/^Ie+i/2 (^) ^-^ ^^^e modified spherical Bessel function of the first kind, JAS9^ . also 



{2£ + 1)J—I,^,/,{k 



2k ' (2£-l) 



4 



see l^j- Hence C is a covariance function with spectrum 

Ji 



Note sinh(K) /k = \/^i\i2 ('^); o.''^^ 



exp (a cos ^9) = ^ + 1) \ —Ii+1/2 (a) (cos i9) (3) 



see USSf , WBMSel . 

Example 6 For k > 0, 7 G [0, vr] ?i;e /iat;e 



00 f 

— (cos 7) = exp (k cos 7) Jo (/« sin 7) 



(.1 



£\2i + V 

see JPBMSBjl, hence exp (KC0S7) Jq (Ksin7) / exp (k) is a correlation function, where Jq de- 
notes Bessel function of the first kind. 

Example 7 We have a covariance function of the form 



C (cos 7) = Iq ( 2k cos 7 — 1 ) Iq ( 2Ky^cos7 + 1 



00 

K 



nPt (cos 7) 

U\? 2i + 1 



see IPBMSSij . Note that Iq (x) is a function see \AS9^ . hence there is no danger of 

complex values. Iq (0) = 1, Iq denotes modified Bessel function of the first kind. 

Example 8 Poisson formula. For a homogeneous isotropic field on we have the spec- 
tral representation 

C{r)= jo (Ar) $ (dX) , 



Jo 

of a covariance function, J Yad83], and the restriction of C (r) on the sphere 'S>2, Co (7) 
C(cos7), is isotropic with spectrum 



^£+l/2(A)l^^(dA). 



(4) 
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Example 9 For a homogeneous isotropic field on M"^^, n > 2, we have the spectral repre- 
sentation for a covariance function 

POO 

C{r)= / i„_3 (Ar) CD (dA) , 
Jo 

and the restriction of C {r) on the sphere S„, C (r) = Co(cos7), is isotropic. The natural 
expansion of C (r), now is in terms of Gegenhauer polynomial C^" ^^^"^ (0037) \Yad83f . 

Example 10 Laplace operator on M^. Consider the homogeneous isotropic field on 
according to the equation 

If m = l, its spectrum (lYadSS^ . n = 3), is 
with covariance of Matern Class 

1 K_i/2 (cr) 



C(r) 



(27r)3/2 2c(cr)~i/2 

1 (cr)^/^ (cr) 



(2vr) 



3/2 2C 



in terms of modified Bessel (Hankel) function, jlASO^ . Now apply ^ for $ (dX) = S (A) dX, 
and we have 



■00 \ 
2 - - 



•^^+1/2 (;^2 + c2)2'^^' 



see WBM86^ . 



3 Non-Gaussian isotropy and the angular spectrum 

In a physical phenomenon the isotropic property is treated as a principle. It means that 
there is no reason to make difference between directions. The corresponding property of a 
stochastic field is that the finite dimensional distributions remain unchanged after rotating 
the space. 
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3.1 Isotropy on sphere 



Definition 11 A stochastic field X {L) on the unit sphere §2 is isotropic (in strict sense) if 
all finite dimensional distributions of {X (L) , L € §2} 0,1"^ invariant under the rotation g for 
every g G SO (3). 

From now on we do not assume Gaussianity. But, for the simplicity, we suppose the 
existence of moments and that those determine the distribution as well. The general form of 
a mean square continuous field is given in terms of spherical harmonics 



where {Z"^; ^ = 0,1,...; m = —£, 1 —£,... ,—1,0,1, ...,£ — is a triangular array, in the 
row £ we have 2£+l elements. Notice that Yq = 1, put Zq = ij., otherwise EZ^^ = 0, therefore 
EX (L) = fi, and the convergence is valid in mean square sense |Leo99j . This development 
follows also from the stochastic Peter- Weyl Theorem, see |MP11) for details. 

In case the m*'^ order cumulants of X are invariant under the rotation g for every g G 
SO (3) then it will be called isotropic in m*^ order. Naturally a strictly isotropic field with 
ru^th Qj-(jgj. moments is isotropic in m^^ order. 

We list here some properties which follows from isotropy, see [BM07] for details, 

1. The distribution of both ReZ™ = ImZ™ and ReZ™/ImZ™ are Cauchy. 

2. Marginal distribution of both ReZ™ and ImZ^ are always symmetric 

3. For fixed £, Zp, m = 0,1, ...,£ — 1,£ are independent iff they are Gaussian. 

Now, let us consider a rotation g G SO (3), it is known that the spherical harmonics 
at the rotated location are given in terms of the Wigner D-matrix, see [TBI Appendix iBj more 




(5) 



e=o m=-e 




A (g) Yr (L) = Di% (g) Yl (L) 



k=-e 



where A{g) denotes the operator according to the rotation g, A{g)Y^ (L) = Y^ ^L). 
Hence the rotated field has the following form 



A (g) XiL) = Y,Y: Zr Yl Di% (g) Yt (L) 



e=o m=-i k=-e 
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The isotropy assumption is equivalent to that the distribution of the variable 



i9) = E (9) zr, (6) 

m,=—i 

is the same as the one of Z^. This statement will be used frequently below. 
In matrix form 

Ze{g)=D^'^ {g)Ze, 



where Zf 



Zj^ ,Zj^ , . . . , Z/^ 



(the dependence on g will be omitted 



unless it is necessary). This equation provides an equation for the cumulant functions 
(log of the characteristic function) as well, in case of isotropy, for each rotation g we have 



(^) = {^D('^ (9)) , 



where lj^ = [a;_^,a;_^+i, . . . ,a;£_i,a;^]. Hence the distribution of Z^ should be rotationally 
invariant on R^^+^. For instance the covariance matrix Cz ((1,^2) = Gov (Z^^ , Zi^) commutes 
with since is uniter and Cz (^1,^2) = D^^^>Cz {ii,£2)D^^^K If ii = £2 = £, the 

only matrix which commutes with D^^) (g), for any g G SO (3) is constant times unit matrix. 

Assume for a moment that Cum2 (Z|;^,Z™) does depend on £i,k,£2,'m, then we show 
below that from the isotropy Cum2(^^,2™) = Se^/,^5^k,ni{—'^)"^ C2 {£1) follows. Hence 
we should restrict ourselves on an uncorrclatcd generating array Z^, with EZ^ZJ^* = 
^e,k^m,n'^£ m- Since o"!^ docs not depend on m we denote it fe, moreover we shall consider 
the field X in the following form 

00 e 

e=o m=-l 

where EZJ^Z"^* = 5£^k^m,nf£- In other words it is seen that the second and higher order 
structure of the generating process inside the same degree £ are hiding, they are not 
identifiable. 

According to the angular momentum of degree i we define the field 

ue (L) = E zrrr (L) . (7) 

m=-e 

We shall be interested in the invariance of the distribution of (L) under rotations as well. 
If the location is fixed then (L) is connected to the distribution of directly, since the 

Condon and Shortley phase convention \f^^^D^Q (g) = Y^* {'d,ip), where the rotation g 
given in Euler coordinates (7, tp) , provides 
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observe that 7 here is arbitrary, see ()32p . Our main interest is the comparison of the finite 
dimensional distributions of the field U£ (L) to the rotated one. In case a rotation carries the 
location L to the North pole A^, U£ simplifies 



Rewrite X into the form 

00 

X(L) = J^u, (L). 

e=o 

We consider a real valued X (L), therefore Z^* = (-1)'" Z^"", since YJ^ (L)* = (-1)"" F"™ (L). 
Moreover if we reflect the location to the center then from {—L) = (—1)^1^™ {L) follows 

00 I 

X{-L) = Y, {-lfzTYr{L). (8) 

l=Q m=-l 

Therefore we have the following 

Remark 12 Since the parity ui{—L) = (— l)^nf (L) is valid, it implies that under assump- 
tion of isotropy 

Gump [ui^ (Li) ,...,ui^ (Lp)) = (_i)^i+^2+-+^p Gump (u^ (Li) ,...,ui^ {Lpj) , 
hence for p > 2 the sum Cp = ii + £2 + • • • + ip must be even. 

3.1.1 Second order Isotropy 

Gonsider the second order cumulants of , defined in ([6]) , 



Gums (Zl , Z^) = ^ti^tl Cum2 [z^ , Zj 



P,g=-<:iv 



Now assume isotropy, Gums [Z^_^,Z^) = Gums and integrate both sides over the 

sphere according to the invariant Haar measure, we have 

Gums [Z£^,Z^^ j - 2^ 2hTl V ^1 ' 

p,g=-£i/2 

p=-£i 

= ^ei,e2^-k,m (-1)™ Cs {£1) 
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where 

^^(^1) = ^ E Cum,(V*,ZfJ. 

Hence Cum2 [Z^_^,Z'^*) = 5i^^i^5k^mft^ i-e- the series is uncorrelated. 
If the series is uncorrelated then 

P,q=-ei/2 

= ^^ui. E Di'^jD(^l{-ir Cum, [Z-^\ZI) 

p,q=-ll 

p=-ei 
= ^ei,e2^-k,m (—1)™ 

since d[^^^^ is unitary, hence Cum2 (^Z^^,ZlJ^*^ = 5i>^^i^5_k,mfh = Cum2 [z^l^^Z^^*^. 

Lemma 13 The field X (L) is isotropic in second order iff the triangular series Z^ is un- 
correlated with variance fi. 

We conclude that a field X {L) with Gaussian i.i.d. Z^ is strictly isotropic. It follows 
from the isotropy that u^^ and u^jWith different degrees £i and £2 are uncorrelated, and 

Cum2 {u, (Li) , u, (L2)) = E Cum2 (Zf ^ , Zp) (L,) Y^^ (L2) 

mi ,m2=— £ 

= E /£^r'"Mii)^r(i2)(-in 

m2=— ^ 

by the addition formula see ([39|) . In particular Var (n^ (L)) = ^^^/^. Instead of the addition 
formula one arrives to the same result if applies the isotropy for Cum2 (u£ {Lf) , (^2))- Ro- 
tate the locations Li and such that becomes the North pole N and at the same 
time Li belongs to the plane xOz. This rotation will be denoted by 5l2Li- We have 
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Cum2 {ui (Li) , ue (L2)) = Cum2 (n^ {gh^L-^Li) , ut (iV)) , since U(, (AT) = ^ and 



Cum2 {ui (Li) , ui (L2)) = Cum2 | ti£ (^LjLiLi) , 



^^"^^ (v ^^^^^^^^ • ^) V 



3.2 Spectrum 

Consider the covariance function C (Li, L2) = E (X (Li) — /x) {X (L2) — /Lt) for an isotropic 
field. Let the rotation gL2Li be the one which takes the location L2 into the North pole 
A^, and Li into the plane xOz. The Euler coordinates of gi^LiLi is the co-latitude an- 
gle -d and 0,since the rotation does not change the angle between Li and L2, such that 
COS1? = Li ■ L2. Under isotropy assumption the joint distribution of X (Li) and X {L2) 
equals to the joint distribution of X (N) and X {gL.^i-^Li), i.e. X [N) and X 0) con- 
tain all pairwise information. In other words the covariance of an isotropic field depends 
on 1? only, Cums (X (Li) , X (L2)) = Gov (X (Li) , X (L2)) = C {gi^L.Li, N), necessarily 
C {Li,L2) = C{Li- L2) = C (cos??). Now 

00 I 
£=0 m=-e 

£=0 m=-e 

and 

00 £ 

X {gL,LM) = E E {gL,LM) ZT, 
£=0 m=-e 

We have Cum2 {Z^, Z]^) = (5^,jfc(5o,n/^, hence 

c (Li, L2) = E H -^Yi {gL,LM) 
1=0 

^=0 
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Similarly to the time series setup the covariance C (Li,L2) is expanded in terms of an or- 
thonormed system ^^^^P^ (cost?) with coefficients In particular the variance Var (X (L)) 
is split up into the superposition of fi's. Hence 5*2 (i) = ft is called spectrum of the field 
X (L) according to the frequency i. 

4 Bispectrum 

If the field X (L) is non-Gaussian then the first characteristic after the second order moments 
to be considered is the third order cumulant, referred to bicovariance or 3-point covariance 
also since it is the third order central moment. The corresponding quantity in frequency 
domain is the bispectrum. Similarly to the spectrum when the covariance (second order 
cumulant) 

Cum2 (X (Li) , X (L2)) = fi^r^' ' 

is split up to the superposition of some orthogonal system of functions and the coefficients are 
called spectrum. We see that the orthogonal and normed system is the Legendre polynomial 
system {^^^^Pt] and the spectrum ^2 {t) = fe corresponds to the frequency £. We put a 
similar question for the higher order angular spectra as well. Namely, under assumption 
of isotropy of the m*'^ order cumulant we consider the series expansion according to an 
orthonormed system and the coefficients will be called m^^ order spectrum. The Wigner 3j 
symbols 

\rni 1712 rrij, J 

will be intensively used from now on, see Appendix [Bl [T21 It depends on the quantum 
numbers ii,i2, ^3 called degrees and orders mi , 1712 , . We list some basic properties called 
selection rules for Wigner 3j symbols 

• 3j-symbol vanishes unless mi + m2 + 7713 = 0, 

• Integer perimeter rule £1 + £2 + £3 is an integer and even if TTii = 7712 = m.3 = 0, 

• Triangular inequality |£i — ^2! < ^3 < ^1 + ^2 must be fulfilled. 

4.1 Isotropy in third order 

The rotation applied to he triangular array shows the necessary and sufficient condition 
for the third order isotropy. 

Lemma 14 The field X (L) is isotropic in third order iff the bicovariance 
Cuma (^Z^^ , Z^' , ^e^^^ '^f triangular series Z^ has the form 
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and 



See Appendix [O for the proof. Note here that Cums (^Z^^ , , Z^^^ = 0, if nii + m2 + 
1713 7^ 0, that is not all the possible bicovariances come into the picture, moreover the cumu- 
lants are depending on the orders 1711,1712,1713 through the Wigner 3j symbols only. In other 
words the third order probabilistic property inside a fixed quantum number £ does not show 
up. The function S3 of the frequencies is an average of the cumulants Cuma (^Z'^^ , Z^^, Z 

by 'probability' amplitudes, see (|26p hence it is called as 'angle average bispectrum', |FSn9) 
we shall turn back to the notion of the bispectrum later in this section. 
Notice first that the quantity 



, S3 (4,^2,^3) 

mi 1712 "23 ' 

is symmetric in £1,^2, ^3, see Appendix [UJO More over B3 (^1,^2; ^3) is symmetric in £1,^2 
as well, since by parity transformation ^ ^1+^2+^3 must be even, the coefficients 



are symmetric as well. 

Prom now on we fix the order of ii,i2, ^3 such that ii < I2 ^ ^; 
Now we turn to the angular momentum field U£. First observe 



ll 12 <^3 

kx k2 k-s 



Cums {ue, (Li) , u,, (L2) , ue, (L3)) = J2 (^1) (^2) (^3) Cum3 (z.^^ 

E ^r(^i)^r(^2)yr(^3)Q ^^3(4,^2,^3) 



mi,m2,m3 



This is invariant both the rotation and ordering of ^1,^2,^3- The 3-product of spherical 
harmonics Yf^, 

.,,,,,3(L„L2,L3)= E (^^ ^)y,-HLi)Y-^L2)Y-HL3) 



is rotational invariant, see Appendix pop. therefore without any loss of generality we apply 
the rotation which takes the location L3 into the North pole and at the same time 
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takes L2 into the zOx- plane 



Cums {ui>^ (Li) {L2) ,ui>.^ (L3)) = Cunis {ui>^ (5^3^3X1) , (5'L3^'2^2) ,-"^3 (^)) 



mi=— ^1 m2=~t2 



mi,m2 ^ 

The third order cumulants of contains an orthonormed system of functions 



2^3 + 1 



47r 

mi;2 



In this way Iii/2/3 i^ij ^2, L^) is connected to the bipolar spherical harmonics, i.e. to the 
irreducible tensor products of the spherical harmonics with different arguments, |VMK88] . 

[Ye, (Li) Ye, (L,)],^„ = (-1)^-^^+- VWTl J] (^^^ ^) (^1) (^2) 

Rewrite Ie^,i2/3 terms of Euler angles 

^^1,^2 A (''9i)V'i)''?2) = Iei,e2,e3 (51.3^2^1! 5L3i2 ^2,^) 



mi,m2 



47r 



According to the usual measure Cl (dLi) {dL2) = sm'i}id'&idipism'd2d'&2d(p2 on £ [OjVr] 
if £ [0, 2-k], we have 

/ Yr 0) y/ (i?, 0) J7 (dL) = 5ej5^^k, 
JS2 

hence 

yy '^ti/2/i^h,h,33^{dLi)^{dL2) = hi,jik2,j2^h,h^ 
§2 

since 



mi ,7712 



(2^3 + 1) Y^r^ ;m =1. 

\ m-i 771,2 / 
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The third order cumulants of ui simphfies 

Cuma {ui^ {Li),ui^ {L2) ,Ui^ (L3)) = B3 {hA2,h)1-h/2/z i:&i,^i:'&2) ■ 
If we denote it by 

Cums (n^j (Li) ,u^2 (^2) ,^£3 (^a)) = C'vi.&A {^1,^1,^2) , 
then we have the inversion formula 

jj C'«,jij2,i3 (i9i;95i,i?2)2:<?i,£2/3 i'&i, ^Pi,'&2) ^ (dLi) n {dL2) 

= ^h,jl^i2,j2^i3,j3B3 (4, -^2, ^3) • 

4.2 Bispectrum of isotropic fields on sphere 

The use of the bispectrum for detecting non-Gaussianity and nonhnearity is well known in 
time series analysis |SRG84] . [Hin82| [TM98] . the similar question has been put and stud- 
ied for CMB analysis see [Mar04j . [AH12j and references therein. The asymptotics of the 
bispectrum has a long history started by |RvN65] . the asymptotic distribution of angular 
bispectrum when the degrees are increasing is considered by |Mar06j . |Mar08j . 

Start with the general case Cuma {X (Li) , X (L2) , X (L3)). The value of the bicovariance 
does not change under the rotation 3x3^2 

Cuma {X (Li) , X [L^) , X (L3)) = Cumg {X {gL.L^Li) , X {gL,L2L2) , X (N)) 

= Gums {X {^1,^1), X {^2,0) ,X{N)), 

the spherical coordinates of the locations gL3L2-^i = ('^ij^'i)) 91^12^2 = ['&2i^) are defined 
by the locations Li, L2, L3 as follows: 5^3X2-^2 ■ N = L2 ■ = cos??2i 9L3L2L1 " X = 
Li ■ L3 = cos-i?!. These angles (-i?!, (^^i, 'i?2) define uniquely, up to rotations, the trian- 
gle given by locations Li, L2, L3. In general the bicovariance is written C3 (i?!, 921, '!?2) = 
Cum3 (X , X (i?25 0) , X (A^)), i.e. it depends on a location L(t?2,0) form the main 

circle {(p2 = 0) and a general location Lq = Lq ipi). 

Cum3(X {Li),X{L2),X{Ls)) 

= E E E Cnms{zr:,Z-\Zl)Y-^{g,,MY,Ti9L3L2L2) 

4/2/3=0 mi=-eim2=-i2 
00 

= Yl Bs{h,£2j3)Te„e2,hi^uVu^2). 
This series expansion of Gum3 {X (Li) ,X (L2) ,X (-L3)) implies the following definition 
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Definition 15 The bispectrum of the isotropic field X (L) is 
and the bicoherence of X (L) is 

■93(4,^2,4) ^ ^3(4,4,4) 

4.2.1 Particular cases 

For any locations Li and L we have 



Cums {ue, (Li) , ue, (L) , u,, (L)) = ^\\^^^ ( q q o) (^i' ^2, 4) ^'^i (cos L • Li) , 
from this readily follows 



h,h,e3=0i=l ^ ^ 

X S3 (4, £2, 4) f"^! (cosL-Li). 

In particular, according to the Condon and Shortley phase convention, we have for any 
location L the third order central moment 

Cuma (n,, (L) , n,, (L) , u,, (L)) = ^11^^ ( q q q') S3 (4,4,4) , 

and the third order cumulant (central moment) Cuma (X (L) , X (L) , X (L)) of X (L) is 
Cum3 {X (L) , X (L) , X (L)) = Cum3 (X {N) ,X{N),X (iV)) 



47r 



t nv^Cj 'J 'j)«3(...,.) 



Note that the summation is taken over those ^1, £2, ^3 when ^1 + ^2 + ^3 = is even and for 
this case the following explicit expression is valid 



h 4 4\ , ,^£/2 / n(^-24)! (^/2)! 
q) ^ ' Y (/: + !)! n('C/2-4)! 
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4.3 Linear field 

First we assume that the rows of the triangle array contain independent variables. In 
other words the angular momentum field is linear. Then for a fixed degree i 

\ t I I J m mj 

now the selection rules apply, see Appendix iBl [T3] 

Cnm, {Zj^\Zf\Zp)=5m^=m5m=o{^^ Q ^^B.iUa)- 

Hence the only nonzero third order cumulant might be Cum^ (^Z^ , Z^ , Z^^ Further the 
bispectrum 

B,{£,e,e) = Y,(J^ I ^)cum3(z,^z|,z|) 

Cumg (ZO, Z", ZO) , 







therefore 



2 

Cums (ZO, ZO, Z°) = (o o) Cums (Z°, Z°, Z°) . 



We conclude from this that from the isotropy and independence follows Cums (Z™, Z™, Z™) = 
0. Moreover if all the members of the triangle array Z™ are independent then 
Cuma {X (Li) ,X (L2) ,X {L^)) = 0. Third order cumulants are vanishes for instance when 
the distribution is Gaussian, or symmetric, etc. 

4.4 Symmetries of the bispectrum 

The cumulants Cum3 (X (Li) , X (L2) , X (L3)) according to different locations ^1,-^2 and L3 
are defined by the spherical triangle with vertices -^1,^2 and L3. The efficiency of computa- 
tions and avoiding the redundancy in statistical procedures require to determine the principal 
domains. The principal domain for the bispectrum according to the frequencies ^1,^2,^3 is 
the following 

1. ^1,^2,^3 is monotone: ^1 < ^2 < ^3, 

2. ii + £2 + ^3 is even, see Remark [T2l 

3. ^1,^2,^3 fulfils the triangular inequality |^i — £2! < ^3 < ^1 + ^2, 

4. Cum3 (z^^^^ , Z;^2 ^ ^rn3\ ^ ^j^^ggg mi + m2 + ms = 0. 
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Similarly the principal domain for the bicovariance according to the locations (Li, L2, L3) 
is {{"&!, <^i, 1^2) G [0, tt] ,(/?2 G [Oi"^] 5 ''^2 G [0, tt]}. We apply the following notation of these 
angles cosi?! = L2 • -L3, cost?2 = Li ■ L3, and ip2 = (j^z is the surface angle at L3. The third 
central angle is given by cosi?3 = Li- L2. The surface angle cps can be calculated for instance 
from the cosine formula 

cos T?3 = cos T?l cos t?2 + siu T?i siu '&2 COS ip2 ■ 

5 Trispectrum 

5.1 Isotropy in fourth order 

In the papers |Hu01| . |KK06) for instance, the notion of trispectrum is used for the coefficients 
of the series expansion of the fourth order moments. Naturally, because of the expression of 
the moment via cumulants (jlSp it contains terms which correspond to the Gaussian part of 
the field. This is one of the reasons that the time series analysis considers the series expansion 
of cumulants instead of moments and reserve the notion of trispectrum for this case only, see 
|Bri65) ■ |MH95j . The cumulants up to third order equal to the central moments, but it is not 
so for higher order, see Appendix |Al Cumulants are very useful, for instance the measuring 
the distance of a distribution from the Gaussian one. 

Lemma 16 The field X (L) is isotropic in fourth order iff the cumulant 
Cum4 of the triangular array has the form 

Cum4 ( Z™^ , ZY"^ , zJ'A = y (-1)™' ( ^1 ^2 \ f i^ h U\ 

xV2i^n{h,e2,h,h\£^) (10) 
See the Appendix iDl for the proof. The function T4 (^1,^2,^3,^4! i^) has the form 

xcum4(4^4^4^4*) (n) 

We notice here, the (jlip shows that the cumulants of Z^ and the function 
Ti (-^1, ^2, ^3, ^4^"'^) define each other uniquely. 

Remark 17 // we assume that T4 (^1,^2,^3,^4! i^) is known for ^1 < ^2 ^ ^3 ^ ^4; then the 
cumulants Cum4 (^Zj^^ , Z^^ , Z^^ , Z^"'^ are given for any {ii ,£2,^3, ^4) • 
Vice versa T4 (^1,^2,^3,^4!^^) can be calculated for any (£1,^2,^3,^4) by <f77]). 
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The sum in (jlip goes for ki,k2, k^, k^ and k^ = ki + k2 = —k^ — k^. Hence the summation 
contains those /cm's when equation /ci + A;2 + ^3 + ^4 = 0, fulfils. Similarly 
Cum4 (^Z^^ , , '^Ti ^ '^TS) ~ ^' ~'~ + "^-3 + "i4 7^ 0. The triangular inequality for 

^i,^2i^^ and ^3,^4,^^ suggests that the 'quadrilateral' with edges (^1,^2,^3,^4) consists of two 
triangles ^1,^2,^^ and £3,£4,£^. The cumulants are invariant of the order of their variables 
therefore this should be true for the other diagonal as well £4, ^i, ^ and £2, ^3, ^- Actually since 
the left hand side of (|lUp is symmetric in ^1,^2,^3,^4, the right hand side does not depend of 
the choice of the diagonal. 

Note that parity transformation ([5]) implies that £1 + ^2 + ^3 + ^4 must be even and we 
turn to the field n^, see [71 

Cum4 (Li) , ui^ (L2) , -Ufa (-L3) , u^^ (L4)) 

4 

m4 



mi:4 j = l 



t2 £i ^ 



^-^ ^ ' Z-^ XL f-j ^ ymi m2 —m \m rns 

The cumulant in the left hand side is symmetric in ^1,^2,^3,^4, and invariant under rotation. 
The function 

4 



(Li,L2,L3,L4) = V2e + 1 E (^^-^ 

h i2 \( h h 

mi 1712 —m}) \m} m^, 1714 

is invariant under the rotation see Lemma [25] Appendix. Therefore we apply the rotation 
5X4^3, which takes the location L4 to the North pole and at the same time takes L3 into the 
plane xOz, 

Cum4 (Li) , (L2) , (L3) , Ui^ (L4)) 

= Cum4 {ui^ {gi^LaLi) , ue^ {9L4,LiL2) , {guLzL^) , ue^ (N)) 



2h + 1 



4tt 



^ ^1 ^2 ^' W ^' ^3 *4 



mi,m2,m3,m'^ j— 1 

we have 

Cum4 {ui-^ (Li) , (L2) , li^a (L3) , n^^ (L4)) 

= (^l,^2,^3,4K^)2:<.^,£2/3/4/l (l9l,l?2,t?3,V5l,V'2) , 
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where the function ^ ^ £ is given by 



mi;3,m}- i=l 



we 


is 


iA 


j \ rn} 







Actually 2^^1/2,^3/4,^^ rotated version of -^^^,^2,^3,^4,^^' corresponds to the tripolar 

spherical harmonics defined as irreducible tensor products of the spherical harmonics with 
different arguments |VMK88j . p. 160. We define the spherical coordinates ("i^i, t?2, "i^s, Vi; V2) 
as follows gi^L^L^ = (f?3,0), gi^L-^Lj = {i9j,ipj), j = 1,2. Note the orthonormality of the 
system according to the measure nfe=i^('^^A;) = Ylk=i^^^''^kd'&kd^k, £ [0, tt], ip^ E 
[0,2^], i.e. 

3 

k=l k=l-A 

^2 

since 

^ ' ^-^ \mi m2 —m l \m U I 

,Tni,Tn2,ms 

here we sum it up by mi, 1712 first then by m^^m^. 

The expression for the Cum4 {X (Li) , X (L2) , X (L3) , X (L4)) is straightforward 

00 

Cum4(X(Li),X(L2),X(L3),X(L4)) = ^ Cum4 K (Li) , . . . , n,, (L4)) 

^1,^2,^3,^4=0 

00 

= X] ^4 {£lj2,i3j4\ i^) 2"^i/2,^3,^4,^l (-^1' . . . ,^4) . 

^i,...,£4=0 £1 

Definition 18 T/ie trispectrum of the isotropic field X (L) is 

Si (^1,^2,^3,41 i^) = n (^1,^2,4,4K') . 

In particular, we have formulae for the marginal fourth order cumulants 



Cum4(X(L0,X(L2),X(L),X(L))= f (0 o) 

X T4 (4,^2,4,4K')2:,i,,2,£i (^i></'i,^2) , 
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and 



Cum4 {X (Li) , X (L) , X (L) , X (L)) 

CXD 

= E 



\j=l 



^1 



also 

Cum4 {X (L) , X (L) , X (L) , X (L)) 



5.2 Linear field 

If the triangular array contains independent variables, i.e. the angular momentum field 
is linear. Then for a fixed degree i 

Cum4 {ZJ^' , Zp , Z^^ , Z^^) = 5m,=m V f ^'i) ( 1 \ 

X (-ir' ^/2i^ + lT4{e,£,e,e\£^), 

from the selection rules follows = — (mi + 7712) = —2m, m^ = ms + m4 = 2m, therefore 
mj = 0. Now for similar reason 

£ £ £^\ f£^ I 
j I 



r4 (£, £, £, £| £^) = V2^i + 1 Cum4 
The only nonzero cumulants are 

cu^,{zlzlzlzi) =Y^^WV\(^^ [ [ ^^\u^,[zlzlzlz^), 



hence for an isotropic linear field U(, Cum4 (Z^^, Z^^, Z^^, Z^") = 0, for and m^. If ad- 
ditionally the series of U(, is independent then Cum4 (X (Li) , X (L2) , (-^3) , X (-L4)) = 
follows. 
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6 Higher order spectra for isotropic fields 



The generalization of the bispectrum and trispectrum is possible for arbitrary higher order 
[MP 10) . [MPllj . If the bispectrum and trispectrum are zero then the field can not be Gaus- 
sian, in case all polyspectra are zero then the isotropic field is necessary Gaussian. First 
we show that the characterization of the isotropy in p^^ order can be given similarly to the 
previous cases. 

Lemma 19 The field 

oo I 
£=0 m=-l 

is isotropic in order (p > 3) iff the cumulant Cum^ (^Z^^ , Z^^ , ■ ■ ■ , Zj^J^ of the triangular 
array has the form 

a=0 ^ 

m} ,...,mP~^ 
P-3 

xllV2FTlSp{ii,...,£p\i\...,iP-^) , 

a=l 

where £^ = £i, i^^'^ = £p, = —ki, k^~'^ = kp, mP = —mi, mP~'^ = rup. The function 
Sp {ii, ...,ep\f,.. . has the form 

sp{ii,...,ip\i\...,i^-^)= (-i)"-^"n(_'I. Izi f-O ^''^ 

ki,...,kp a=0 

k'^,...,kP-^ 

P-3 

X J] VWTlCump [zl\Z^^,. . . , Zll) . 

a=l 

See the Appendix [E] for the proof. Consider a convex polygon with vertices Ai-p-i, and 
edges The diagonals denoted by j = 1, 2, . . . ,p — 3, starting from the vertex Ai 

divide this polygon into p — 2 triangles. The first triangle has sides (angular mometum) 
£2 and i^, the next one has sides £3 and i"^, the general one is £"', £a+2 and i^^^, finally 
the last one £p~^, ip-i and £p. For each a the sides of the triangle [£"",£ a+2, £"'~^^) should 
fulfil the triangle inequality |r -^"+^| < 4+2 < + ^"+^ The coefficients in 1^ will 
differ from zero if orders —m'^,ma+2,^'^^^, fulfil the assumption —m'^ + ma+2 + m'^^^ = 0, 
for all a. This implies mi + 7712 = —m^, m^ + m^ = m^, . . . , mp-i + mp = m^~^. Let us 
plug in consecutively m° and we shall arrive to the result mi + m2 + • • • + mp = 0. Hence 

Gump (^Z^^ , Z^^ , • • • , Z^^^ = 0, unless mi + m2 + • • • + mp = 0. 
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Remark 20 The cumulant Currip yZ^^ , , . . . ,Z^^j is invariant under the order of the 
quantum numbers li,i2, ■ ■ ■ ,^p, hence the right hand side of is invariant as well. The 
result is that the function Sp . . . , £p\ £^ , . . . given on values £i < £2 ^ ' ' ' ^ £p will 

determine all cumulants Cuirip (^Z^^ , Z^'^ , . . . , Zj^''^ by U^) . 

Repeat the representation of the field 

00 

where 

e 

ne{L)= Yr{L)zr. 

m=—l 

We are interested in the invariance of the distribution of ui (L) under rotations. If the location 
is fixed then u^{L) is connected to the distribution of Z'l, since the Condon and Shortley 
phase convention (f33|) provides 

m=-i 

Our main interest is the comparison of the finite dimensional distributions of the process 
ui (L) to the rotated one in case the rotation carries one location to the North pole N , since 
U£ simplifies to 

»'<"> = \/^< 

We have 

Cumj, (n,, (Li) , u,, (L2) ,...,ut, (Lp)) = E H ^ ^^^^ ^um,, [z^^\ Z^^\ Z™-) 

mi,...,mp j=l 

mi,...,mp j=l i^,...,£P^^ a=0 

,...,mP~^ 

p-3 

xl[V2FTlSp{£u...,£p\£\...,£P-'^) 

a=l 

= E hi.,p/^-p-3 (Li, . . . , Lp) Sp {£1, . . . ,£p\£^, . . . ,£P~^) . 

£i,...,£p-3 
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The following p— product of spherical harmonics 

_ p 
hr,.../p/\.../p-s (Li, . . . , Lp) = Yl ^i^' i^j) 

mi,...,mp j—i 
•m} ,...,mP~^ 

11 \ _j^a fn^,2 m"'^^ I J-J- 

a=0 a=l 

is rotation invariant, see (f30|) . hence we can apply the rotation gipLp-i such that gipLp-iLp = 
N, and gipLp^iLp-i belongs into the plane xOz, we have 



^£i,...,£p,£i,...,£p-3 (^1, ■■■,Lp) = Y n {9LpLp.,Lj) 

mi,...,mp j—i 
,...,mP~^ 

^ ' y-m^ ma+2 m°'^^J J-JL y-mP mp_i Oy ^ ^ 

According to the usual measure 11^=1 ^ (dLk) = Y[^=\ si^^kd'&kd^k, '&k £ [0, vr], (/j^ € [0, 27r], 
the system of functions l£-^^,,,/p/i^,,,/p-'i forms an orthogonal system. 
Consider now 

mi,...,mp a=o ^ ^ ^ 

X n Cump , . . . , 

a=l 

where remember £^ = ii, and l^"'^ = £p. The left size is symmetric in ii, . . . ,£p hence the right 
size must be symmetric as well. One might fix an order for the entries of ^1, . . . , to get a 
unique representation for the cumulant. We consider a monotone ordering ii < £2 ^ ■ ■ ■ ^ £p 
and refer to it as canonical representation. 

Note that parity transformation ^ implies that ^1 + £2 + ^3 + • • • + must be even. 

00 

Cump {X {Li),X (L2) , X (L3) ,...,X (Lp)) = ^ Cunip [ue^ (Li) , ui^ (L2) ,---,Ui^ (Lp)) 

h:p=0 

00 

^ X] 'S'p (^1, • • • ,-^pK"^, • • • ,-^^~^) -^£i,.../p/i,.../p-3 • • • , ip) • 

ei:p=0 

Definition 21 The p^^ order polyspectrum of the isotropic field X (L) is 
Sp (-^1) • • • ) £p\ £ 1 ■ ■ ■ iP' — Sp (^1, . . . , £p\ d- 1 . . . . 
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6.1 Linear field 

Let us consider the particular case when ZJ^ are independent if i is fixed. We have 



... P'~'^ "=0 
m-'-,...,mP— ^ 

P-3 
a=\ 

We have seen that for p = 3, 4, all mj = follows. Let us see the case p = 5, 

\m\.2 vn} J \y—vn} 7713 w? j \—rr? ni^J 

then from the selection rules follows that m} = — {mi + 777-2) = —2m, m? = m^ — m^ = —3m, 
m^ = 777-4 + w-s = 2777, heuce 777 = 0. In general it is easy to see that 777^^ = — (A; + 1) 777, for 
A; = 1, 2, . . . ,p — 3, and at the same time 777^"^ = 2777, hence 777 = 0, and m^ = as well. 
Consider the polyspectrum 

ki,...,kp a=0 ^ ^ a=l 

k^,...,kP-^ 

from the independence follows that ij = I and kj = k, for j = 1,2, ... ,p. Hence similar 
argument to the previous one leads to the result: kj = 0, and A;-' = as well. We have 

S, {i\ i'--^-') = n (^J I ^ Q^') n V2FTTCum, . . . , Z°) . 

a=0 ^ ^ a=l 

Now, by the Lemma [19] we get 
Cum,(zO,Z°,...,Z°)= Y n(o n^^2FTlCum,(Z°,Z°,...,Z°). 

Hence Cump (^Z'^^ , Z^'^, . . . , Z^^^ = 0, for all kj. Now we have a general conclusion because 
the only case when all cumulants vanish except the second order one is the Gaussian. Once 
the rows of {.^™} are Gaussian then all the entries of {Z^^} are independent. Lideed the 
isotropy implies that all the entries of {Z^} are uncorraleted now they are Gaussian hence 
they are independent. 

Lemma 22 // the isotropic field X (L) is linear, i.e. inside the rows of the generating array 
Z"^ all random variables are independent, then the whole array contains independent entries 
and X (L) is Gaussian. 
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7 Construction of isotropic field 

The stochastic isotropic field 



e=o m=-e 



on sphere has very spacial angular spectra, in particular the form of the cumulants of the 
triangular array {Z^^} is restricted, see (|12p . If one starts with a non-Gaussian continuos in 
mean square X (L) then the triangular array {Z^} naturally given by the inversion formula 



Z^= [ X{L)Y^* {L)dL. 

7§2 



Now, we address the reverse question of construction of a triangular array {Z^} with the 
desired properties of their cumulants. Let us start with triangular array assume it is 



Zjf ,Zjf , . . . ,Zjf 



uncorrelated and all moments exist. Consider the vectors Zg 

0, 1, 2, . . ., according to the rows of The finite dimensional distribution is characterized 

by its cumulant function (logarithm of the characteristic function) 

{uji\ £ = 0,1,2, ...) = In Eexp ^"^ujjz^^ , 

such that only finite many coordinates of the variables {uje\i = 0, 1,2,.. .) are different from 
zero. Consider the transformed series Ze = D^^^Zi, where is the Wi gner matrix of 

rotations D^^^ 



k,m 



, and define a triangular array {Z^^} through the cumulant 

k,m=—£ 



function 

^ziuJe\i = 0,1,2,.. 



[ In E exp ( i y loJd^^'^ Ze ] dg 
JsOi3) V V / 



where again only finite many coordinates of the variables {u£\ i = 0,1,2, .. .) are different from 
zero and dg = sin'dd'ddipd'y / Sir"^ , is the Haar measure with unite mass. This new triangular 

array {Z^} will be called Wigner D-transform of .The third order cumulants of 

for instance 



Cum3«,Z,^^4^) = / E <iAl%Al%.dgCnrn,{z-^,Z-^Zl 

J SO{3) mi,m2,m3 

^(h (2 e-A y (h (2 ^ A Cum- J Zr\Zp,Zp 

(4,^2,^3), 



mi,m2,ni3 

h t2 

ki k2 fca. 
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which fulfils ([H). The function (^1,^2,^3) is given by 

mi,m2,ms ^ ^ 

and also 

The conclusion is that a subset of cumulants Cunis {z^^ , , ^^^^^ > i-e- + m2 + ms = 0, 
is used in the construction and the superposition with 'probability' amplitudes is applied. In 
general we also have 

a„„,(z,7..z™,...,z-)= E (-i)=^""'"°n(I. 

£l^,,,^£P-3 a=0 ^ ' 

,...,mP~^ 

P-3 
a=l 

where 

s,(A,K'-=)= E (-')^"'""r[(_i; „t: 

mi,...,mp a=0 ^ ^ 

m-'^,...,m''~'^ 

P~3 

X n V2FTlCump 4^ . . . , Z,';) . 
0=1 

A Cumulants 

A.l Basic properties 

Cumulants (also called invariants) are very important quantities for the characterization of 
a random series, see [RST06] for some details. Moments and cumulants are equvivalent, 
since the moments can be defined as the coefficients in the series expansion of the charac- 
teristic function similarly cumulants are the coefficients in the series expansion of the log- 
characteristic function. The main importance of the cumulants is that it does not contain 
the Gaussian part of the moments. Let X G i?" denote a random vector. 

A. 1.1 Symmetry 

Cum(X) = Cum(Xp(i^„)), X G M", 

where p(l : n) = (p (1) , p (2) , . . . , p (n)) , p G and denotes the set of all permutations 
of the integers (1 : n). 
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A. 1.2 Multilinearity 

For any constants c\-_2 = (ci,C2) 

Cum(ciyi + C2Y2,X) = ci Cum(yi,X) + C2 Cum(y2,Z). 

A. 1.3 Independence 

If X G i?" is independent of y G BJ^ where n, m > then 

Cum(X,y) = 0. 

In particular if n = m then for independent X and y 

Cum{X + y) = Cum(X) + Cum(y). 

This formula is the additive version of the formula of the moment of the product of indepen- 
dent variables. 

A. 1.4 Gaussianity 

The random vector X G is Gaussian if and only if for all vector with elements from 

(1 : n) 

Cum(X,^^^^^) = 0, m>2. 

A. 1.5 Expression of the cumulant via moments 

The first order cumulants is the moment 

Cum(X) = EX, 

the second and third order cumulants are the second and third order central moments 

Cum(Xi,X2) = Cov(Xi,X2) (14) 

3 

Cum(Xi,X2,X3) = E[](Xi-EX,). (15) 

i=l 

The general formula is 

n m 

Cum(X) = ^(-l)-i(m - 1)! n E X^n\ X G M", (16) 

m=l UGVn j=l 

where the second summation is taken over all possible partitions U, where [uj^i-n]j^i.^ is an 
indicator of a partition C with m subsets of the set (1 : n) = (1, 2, . . . , n), uj^k = 1 if A; G £, 
otherwise 0. The double sum in (|16p is over all partitions of £ G Vn, where Vn is the set of 
all partitions of the numbers (1 : n). 
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Example 23 First note that if n = 3, then the m*^ order partitions C & V3 with the corre- 
sponding indicators and permutations 0/ (1 : 3) are 

1. form = l; C = {(1 : 3)} 

C ^ ni,i;3 = [1, 1,1], E X'^^,'-' = 

2. for m = 2; A = {(1), (2, 3)} , £2 = {(1, 3), (2)} , £3 = {(1, 2), (3)} , 



£1 o 



£2 ^ 



£3 o 



■"1,1:3 




" 1 





" 


■U2,l:3. 







1 


1 


^^1,1:3 




" 1 





1 " 


''^2,1:3 







1 





^^1,1:3 




" 1 


1 


" 


.""2,1:3. 










1 



J]EX,"f3^^^ = EXiEX2X3, 

2 

nEX"f^ = EX2EXiX3, 



2 



= {(1 


), (2) 


(3)} 










?^l,l:3 




" 1 





" 


£ ^ 


^2,1:3 







1 







.^2,1:3. 










1 



n E x";^'' = EX3EX1X2, 



J]EX,"f^^ = EXiEX2EX3, 



Now we apply Uu\) : 



Cum(Xi,X2,X3) = EX1X2XS - EX1EX2X3 - EX2EX1X3 

- EX3EX1X2 + 2EX1EX2EX3 

3 

= El[{Xi-EXi), 

i=l 

The first three cumulants, see also equal the central moments but this is not true 

for higher order cumulants. One might easily check this for the case of cumulants of 
order four. Indeed, X 

Cum (X) = Cum {X - EX) 



El[{Xi- EXi) - Gov (Xi, X2) Gov (X3, X4) 



i=l 



Gov (Xi, X3) Gov {X2, Xi) - Gov (Xi, X4) Gov (X2, X3) 
4 

/eJ](x,-ex,), 

1=1 



(17) 



unless all covariances are zero. 
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A. 1.6 Expression of the moment via cumulants 

We consider the moment Ej|"^^Xj as the general case because the moment EY^.^ can be 
put into the form E]^"^^ Xi, where we define n = and 





Now 

n 

^11^*= E n Cum(Xb^,), (18) 
i=i ceVnhjec 

where the summation is over all partitions C = {bi, b2, . . . , bfc} of (1 : n). 
An example of the use of the formula (jlSp can be seen when n = 3. 

Example 24 Each partition of {1 : 3) are listed in the previus example therefore if X gW^, 

EA1A2A3 = Cum(Ai:3) + Cum(Ai) Cum(A2, A3) 

+ Cum(A2) Cum(Ai, A3) + Cum(A3) Cum(Ai, A2) 
+ Cum(Ai) Cum(A2) Cum(A3). 

Now in particular 

E A^ = Cum(A, A, A) + 3 Cum(A) Cum(A, A) + Cum(A)^ 



B Spherical Harmonics 

The following notation and results are used in this mauscript. 

1. M'^ is Euclidean space of dimension 3. 

2. §2^ sphere with radius 1 in M^, colatitude (angular) coordinate i? G [0, tt], longitude 
coordinate if £ [0, 27r], North pole N: ■!? = 0, (/? is indeterminate, (latitude coordi- 
nate 7r/2 — € [— 7r/2, 7r/2], North pole: = 7r/2,) The spherical coordinates on §2- 
(sin cos ip, sin i} sin (p, cos i9) 

3. SO (3) denotes the 3D (special orthogonal) rotation group 

4. Addition Theorem. [EMUTSlj vol. 2,7.15 (30), pp.116 Let Ui and U2 be two vectors 
with angle and lengths ri and r2 respectively, and the Euclidean distance is denoted 
by P = — U2\\ = + r2 — 2rir2 cos 7 , then we have 

00 

p-'Ju {p) = P^y'^T (u) J] (^ + zv) (cos 7) Ji+u in) Ji+u (r-2) , (19) 

e=o 

where Ju denotes the Bessel function of the first kind, Gegenbauer polynomial 

(c,"^ = p,). 
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5. Lapalce-Beltrami operator Ab = Vb) 



A; 



A/ 



A; 



sm 



d 



1 92 



sin-!?—— + 



d 



d 



1 9- 



i2 n 



d'd J sin (^(/j^ 



92 CDS'!? 9 1 9^ 

+ ^-7 7777 + 



O'f?^ sin I? d'd sin^ i? ^(/^^ ' 
( |Apo07| , vol. 2 pp. 293). 

6. V is the central difference operator (yf + vi) ■^j,k = -^j+i,k + -'^i-i.fc + ^j,k+i + 



7. Spherical function / of order i. 

Ab/ 



1 d 



sin I? 



5/ 
5^9 



+ 



sm 



sin 5'!? 
= -iii + l)f, 

|VMK88] . Regular spherical harmonics Yg, homogeneous polynomial of degree £ and 

v^y^ = 0. 

8. Normalized Legendre polynomial Pi (Rodrigues' formula) 



21 + 1 1 dUx'^-l) 



2 2^ei dx^ 
(^Piix)] dx = l 



xe [-1,1] 



Standardized Legendre polynomial Pq (x) = 1, 

, , Id' (x' - 1)' , , 

(1) = lf |EMUT8lj . vol. 2, pp.180) it is orthogonal and 



[Piix)Ydx 



21 + 1 

27r f'TT 



[Pi {cos i})fn{dL) 



Jo 



[Pi (cos "i?)]^ sin "ddddip 

2tt I [Pe{cos'd)fsmM'd 
Jo 

2tt j [Pi{x)fdx 



4lT 

2£ + l 



(20) 
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recurrence 

{£ + 1) P,+i {x) = {2i + 1) xPe (x) - (x) 
generating functions f |EM()T81j . vol. 2, pp.183) 

oo 

Y,Peix)z^= {l-2xz + zY^^\ xG (-1,1), \z\ < 1 (21) 

£=0 

oo ^ 

^ -P, (cos 7) = e^-^ Vo (zsinT?) (22) 
e=o 

9. Pp is the associated normalized Legendre function of the first kind (Gegenbauer 
polynomial at particular indices) of degree £ and order m 

recurrence formula is 

(£ - m + 1) PP^ (x) = {21 + 1) xPp (x) - (£ + m) P^^ (x) , 
in particular P™ (1) = 6m,o, -P^ (x) = (x), 

10. Funk-Hecke formula ( |M166] ) Suppose G is continuous on [—1, 1], then for any spher- 
ical harmonic (L) 

/ G{Li-L)Yi{L)n{dL) = cYi{Li), 

c = 2Trj G{x)Pe{x)dx, 

where Q (dL) = sin^d'&dip is Lebesque element of surface area on §2; Li ■ L2 = cos'&. 
In particular 

/ GiL,-L)Yr{L)n{dL)=geYr{Li), (23) 

gi = 2tt J G (x) Pi (x) (ix 
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see |Yad83| pp.72, for a general dimension d instead of Pi one has the Gegenbauer 

polynomial ^ (^1^^ — ^ij ■ 
Also 

/ G (Li • L) Pe {L2 • L) n (dL) = f,P, (Li • L2) 

fi = 27r j G{x)P,{x)dx 
here fi = {LT)^ G is the Legendre transform of G. Hence 

G (cos = V / G (x) (x) dxP<. (cos i?) 

— ^ 2 1 _^ 

£=0 ^ 

In particular Funk-Hecke formula when N denotes the north pole 
/ [Pi{N ■L)fn{dL) = cPdN -N) 

c = 2tt j ^[Peix)fdx, Piil) = l 



[P, {N ■L)Yn (dL) 



2£ + l 



11. Orthonormal spherical harmonics with complex values Yp {^,ip), i = 0,1,2, .. 
m = —i, + . . . — 1, 0, 1, . . . , ^ — 1, £ of degree i and order m (rank £ and projection 

m) 



Yr (^, ^) = i-ir y ^fr^jl^r (cos ^) e'-^ G [O, 2vr] , ^ G [0, ^] (24) 

Yri^,ipr = {-irYf"^{i),^). 



In particular y,o (i9, = ^^P^ (cos^9), Y,^ 99) = ^J, 



Yr{N)=5rn,o\l^^.. 



we have (E 



/'27r /"TT 

/ / |yr(^,<^)l^sinMi?d(^ = 1 
JO Jo 



note that Y^™" is normalized fully, some authors do not apply l/v47r in the definition 
of Y^™", also for a sphere with radius P spherical harmonics are normalized additionally 
Ypi^,^) /R. 
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12. Wigner 3j-symbols, |Lou06| notation 



m-i-.s] I mi 1712 rris 



Clebsch-Gordan coefficients, 



^ \ ( ^ ^^3+m3+2^l 



^1:3 A ^ t_Y\h+i2+t3 ( ^2 ^1 ^3 

mi:3y \m2 mi m-s 

h h h " 

7712 ^3 TTT-i 

^1:3 \ ^ /_-|^N<?i+fe+^3 ^1 ^2 ^3 

mi-sj V"'"^! ""1-2 — "l3 



Orthogonality relations 



mi;3 ^ ^ 



-'-^ \mi-2 ml \mi-2 kJ 

r?ii;2 ^ ■ / \ • / 



i,k 



Kl «2 , p 

' ' 



if jO =ii + l2 + l^ is odd, otherwise see |Edm57j . (3.7.17) 



■1 12 i3\_, ,,£/2 /n(^-2^.)! (^/2)! 



o; ^ ^ Y (£ + 1)! n('C/2-£j)!' 



if £ is even let i = 2k, say then 



_r_n3fc / (^0' (3fc)! 

o; ^ ' y (6fe + l)! (fe!)3- 
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13. Selection rules: 

• 3j-symbol vanishes unless nii + m2 + ms = 0. 

• Integer perimeter rule £i + £2 + ^3 is an integer (an even if mi = m2 = = 0) 

• triangular inequality |^i — ^2! < ^3 < ^1 + ^2 

14. Gaunt series [ Edm57) . (4.6.5) 



Mo o)imi.2 m)^^ 



|VMK88] . pp. 144, products of three and more 

15. Rayleigh plane wave expansion in 3D: k = k/\k\,'x = x/\x\, r = \x\, k = \k\, 

00 e 

^ik-x ^ ^ikr cosH = 47r ^ ^ i^je (kr) Yf" (tj * Yf" (x) 

i=0 m=-i 

where 

is the Spherical Bessel function of the first kind, also 

00 

e=-oo 

00 
e=-oo 

see |Lou06j . 

16. 2-D Dirac delta on sphere 

6 {L,,L2) = Yr (Li) YP (L2) = ^ E + 1) " ^2) , 



1-D Dirac delta 

^ (x — 11) = ^ 



i 



see |Lou06] . 
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17. Rotational invariant functions. The function 

I 

21 +~ 

e 



k=-e 



A-TT 

= ^7-lE^/(^i)^^'*(^2) (29) 
k=-e 

is rotational invariant. This expression is valid when we apply rotation g^-^, where 

= Pi (cosLi • L2) . 
It follows from the addition formula (1391) as well. The function 



/,,,,,3 {L,,L„ L3) = ^^^J^^l E (^';3) (^1) (^2) (^3) (30) 

of three location is rotational invariant, see |Lou06j . pp.14. We repeat the notation 
^^■^ 1 = 1^^ I . In particular we have 

"1.1:3/ \mi 1712 m^J 

18. Wigner D-matrix A (g) (L) = {g'^L) 

e 

A (g) Vr (L) = Di'l (g) (L) (31) 



k=-e 



if £ is fixed D^f^ (g) is unitary 



e 

E Dlil*k (a) = ^m„m„ 

k=-i 

:=-e 
e 



k=-£ 
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m,k 



for fixed rota- 



. We introduce the notation D^^^ = D 

tion g. Thus D^^^ denotes a unitary matrix of order 2£+l, and it follows -D^^^ [D^^^] 
L>(^)L»(^)*, unimodular detL>(^) = 1. For the definition of dI^\ (g) 



Dm,k 7) = exp i-imip) exp (-^^7) 



(32) 



where d!"^f^ is the Wigner (small) d-matrix it is real and 



-i)'~"^dS,_(^) 



(^) 

-m,—k 



(^) 



see [VMK88] . pp79 for details 

The notation * defined as the transpose and conjugate for a matrix and just conjugate 
for a scalar, for instance 



D 



m,k 



(-7, -i?, -(/?) 



note here that if the Euler angles {ip, "i?, 7) corresponds to g then the one (—7, — "!?, —^p) 
corresponds to g~^, hence 

dZ is) = Dil {9-') 

also 

dZ = DZn+i,i-k+i (9) 
Let g the successive application of gi and 52 then 

DZi9)=i:DZi9^)DfUg,). 
j=-i 
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Condon and Shortley phase convention, |Edm57| . (4.3.3) 

vr (•'.»') = (7, «, V) (33) 

= (-iroS (7. 

where 7 is arbitrary angle. This form is referred to as passive convention as well. It 
also follows 

Yf"" (i?, if) = (-1)™ e-*2™^yr i^, ^) ■ 



|Wor) ■ (52) in terms of Wigner 3j-symbols, |Edm57| (4.3.2) Singly coupled form, Clebsch- 
Gordan series: 



£,m,k 

^il ,mi ,m2 ^ii ,fci -,£2^ ^m,k ' 



£,m,k 

£l+i2 



^m^M^m2,k2- 1^ \^mi;2 m)\k^..2 k)^ ^> ^^'^ + ^}^ 

£=\£i~£2\ 

—m = mi + 7712, —k = ki + k2, j Edm57j . (4.3.4). Double coupled form 

^ lmi:2 m) ™2,fc2 fcy 2^ + 1 '"'^ 



mi+m2,fci+fc2 



mi,r?i2,A;i,fe2 

|Edm57| . (4.3.3) 

3) /^l:3\_M:3 



y f 

/ ^ mij/ci m2,fe2 msjfea \ 



mi ,m2 ,m3 
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The integrals dg = sin 'dd'ddipd'^ / Svr 



2 



^m.k^g = ^efl^mflSkfl (36) 



m.k 

50{3) 



lsO{3) ^Slfci^Sk^^ = ^h,h^m^,m2h^,k2^l^--[ (37) 



Ml/2,h _ / r^(ei) T^{h) , _ / «1:3 \ / ^1:3 \ .ooN 



J SO{3) "^'^'^'^ rn2,k2 mz,kz » -)- 1 ^l,mi;^2,m2 ^1, 

where the Haar measure is dg = sin ^d'ddipd^i I Stt^ . |VMK88j . 
19. Scalar product for f,g 



2-77 /'TT 



(/,5) = / / f{^,^)g{^,v)sinMM^ 
Jo Jo 

20. Addition formula [GROO]GR 8.814. |EM()T8T] vol. 2, 11.4(8) pp. 236, [Ml66] . 

.P,(l) = l, COS^?=((l9i,(/Pi)-(l?2,¥'2)); 

Of _|_ 1 

y-(^i,c^i)y-(^2,^2) = ^P,(cos^) (39) 

m=—£ 

[Yad83| For general dimension d |Yad83j p. 72 provides an addition formula in terms of 
Gegenbauer polynomials Cp, [C]^'^ = Pi 



m3\ I «1:3 ] ^ ( ^1 ^2 h 
"1-1:3/ 1712 1^3 



C Proofs, Bispectrum 

We use the following notations Z^^^^' = 
Proof of Lemma. [2] 

Cumg (Z^'A = (^'■■' ] (£l:3) , 



mi;3 
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Assume then 

mi, 1712,1713 

= cum3 

= (^Z^^ , Z^^ , Z^i^^ see Edmonds (4.3.3) !![Edm57], [CMOlj . hence the assumption of 
isotropy fulfils. If ^ is not assumed then under the assumption of isotropy we have 
Cums (z'^^;A = Cuma (z'^':^Y and integrate both sides of gOD, see dM]), 




where 

«3W:») = E (,'„';:) Cum, (zs-)- 

Hence ([9|) is a necessary and sufficient assumption for the third order isotropy. In this case 
the bispectrum is a linear combination of the cumulants fo the angular projections by the 
probability amplitude of coupling three angular momenta ■ 
Proof. Consider the expression 




Any permutation of £i:3 implies the same sign of the Wigner coefficients hence it does not 
change the value. ■ 

Proof. 3-product of spherical harmonics is rotation invariant, indeed 
see ([35]). ■ 
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D Proofs, Trispectrum 

Repeat the notations £1-4 = (^1, ^2, ^3, 4), 

Cum4 {X {Li),X (L2) , X (L3) , X (L4)) = Cum4 {X {Li.,4)) , T4 {h , 4, ^3, 4| i^) = ^4 {£i-a\ e^) ■ 
Lemma 25 

V L>(^l) Z)(^2) r)(^3) ^m ( h:2 ^ \( ^ ^3:4 

A:2,m2 fcs.ma fc4,m4 \m1.2 -mj\m m^^-i 

mi:4 ^ ^ ^ 

/'«1:2 « \ ( tSA 
^ ' yh..2 -k \k3A k 



Proof. Indeed, from 



and from mi + 777-2 + ^^3 + "^-4 = 0, follows 



V-r){^i) 75(^2) ^(^3) r)(^4) ._-,Nm/^4:2 ^ ^ ^ ^3:4 A 
fci,mi fca.ma fcs.ms ^4,^4'^ ^ Wl-2 -m \m 7713-4, 

x(-l)-f^3:4 ^'W4:4 ^'W4:2 M ^ ^3:4 \ 
^7713:4 ?77^y V^3:4 k"^ J ^7771:2 — 777^ ^777 m^-^J 

= E Y.{''-' '\){' '^^0(2^^ + 1) 

XI Z^^fci,mi^A:2,m2^fc2,™2i ij ^^^^^^ -fc^y) ^7771:2 "Tnyl 



'^3:4 i 


92 ■ 




,^3:4 - 


k\ 


)(. 








^ I 7771:2 




-777^ 



m,-m ki,Tni fc2,m2 -k^,m? \mr>.n —m.^ I ^ ' 

4:2 ^ /^^3:4 ^ 

fcl:2 ^ ^3:4 



,^4:4 
^3:4 



Now use 



see (123). 



^ V "1-1-2 —777/ V 7771-2 "1- / 

mi:2 ^ ' / \ - / 
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Lemma 26 



J^iil) Jjik) J^ii4) r 

50(3) ''"^ fc4 ,m4 



= j;(2f+i)(-i)-ff- n 

Proof. Here the we apply the formula p4p for the single coupled D-matrices then the integral 
of triple product of D-matrices [38] and get 

^(^1) (.4) ^ 

fel,mi k2,m2 kz,mz k4,m4 ^ 



SO{3) 

/ ^ V^3;4 /i^V V"^3:4 JsO{3) 



-^fci ,mi ^fea ,m2 fci , -ml "y 



(2^1 + 1) (-1)^ 



^1:2 TT^"*^/ V"^3:4 — m"*^ 



Lemma 27 // + £2 + + ^4 is even then 
^(2£ + l)(-l) 



fc-m (h:2 i \ n hA\ (il:2 M / ^ ^3:4 



is symmetric. 

Proof, of Lemma [271 Start with 



il:2 t \ (i i3A 
ki:2 -k) [k ks-A 



^(24 + l)(-l) 
4 



T;+£o-£-ki-k3 ( h h 4 



^2 ^3 —ko J \ki kQ k. 



1 ^0 ^4 \ ^1 ^2 



4/ l<i3 <i4 



where T, = ii + £2 + is + ii, ■ Now we evaluate the product 

-e-k,-k, fi2 is ih \ fh ^0 h\ f ii i2 



J](24 + l)( 1)^+^0 -fciyV^l ^0 ^4/ L.3 .4 .0 



V(2^g + l)(-lf+^o-^-— % ^01? . 

^0 



^1 ^2 * 

<2 
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The values 



{2i + 1) {2el + 1) 



4 4 il 



form a real orthogonal matrix, see |Edm57] (6.2.10), rows and columns being labelled by i 
and Iq. Hence if we sum it up by the result is 



^4 
1714 



q.e.d. ■ 

Proof of Lemma I16L We have 



Cum4 {Z^:;:) = E Cum, (Z,^ 



7m,l:4 
mi:4 



Under the assumption of isotropy Cum, ^Z'^^'^^ = Cum, ^Z'^^'-^^ , now integrate both sides 
by the Haar measure and get 



Note that each term according to summation m-i;, is symmetric in -^i;,, see Lemma [271 
above. Now define 



,mi:4 

with this notation we have the cumulant in the form 
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If instead of isotropy ()10p is assumed then 

mi:4 

n(^2) n(^3) r,(^4) 

A;4,m4 

mi;4 



X 



/ ^ fei,mi fe2,m2 K3,m3 k 
mi;4 



= J^T4(£l:4K)(2^ + l) 

xV-^(^i) r^^m ( h:2 ^\(^ ^3:4 

/ ^ fci,mi A:2,»Tt2 fc3,m3 ^4,014 V J I ?Tii-2 —m J \m TTl^-^ 

mi;4 ^ ■ ' ^ 

the Lemma [251 can be apphed and get 



Cum4 ( 



Particular cases 



Cum4 {u,, (L) , ^z,, (L) , n,3 (L) , u,, (L)) = r4 (^i, ^2, ^3, ^4| 11 (^i 



■J 

ml /I mi Aj=l 

I mi m2 —m} J \ m} ms J ' 



Cum4 (Li) , Ui^ (L2) , ti^3 (L) , (L)) 

= Cum4 (n<?i {gLL2Li),Ui^ {gLL2L2) ,ue.^ (N) ,Ui^ (N)) 



^ i=3:4 



£1 



3 <^4 





\/2^i + 1 



5] y,7(9LL2ii)i^r(5^^2^2) 



mi ,m2 



^1:2 

mi:2 



47r 



^ i=3:4 



£1 



^3 4 





T4(£l,^2,4,4K')%,£2/i (^?2,V'2,^l), 
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Cum4 {ue-, (Li) , (L) , ue^ (L) , ue^ (L)) = Cuir^ {ue^ {gLLiLi) , (iV) , ue^ (iV) , (iV)) 



1^- 



—2:4 ^1 \ / \ / 



^ 



n ^^E^4(<"<2,«3,<4l<')\/5rTTp„(L.Lo(o' '^)( 

— 1 .4 «1 \ / \ 



j=l:4 ^1 



M 



Cum4 (n^j (L) , ue^ (L) , tx^g (L) , ue^ (L)) = Cum4 (u^^ (AT) , (AT) , ue^ (AT) , u^^ (AT)) 



i=l:4 ^1 



E Proofs, Polyspectrum 

Notation: gi^'" = ^f''!''-'!^'' 

Lemma 28 Define f = h, £P-^ = ip, A;° = -ki, kP'^ = kp, m° = -mi, ■wP-'^ = rup 

_ f fTl)(^<^) da- , -^^S(ml^^-3-fcl:POO-3)^' / T 4+2 
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Proof. The fourth order key formula is emphasized here for understanding the induction 

Gt' = f Di'^^ D't^ ) ) D't^ dg 

^fel:5 ;mi:5 / „^ s"^! ^2 ,^2 Ks ,7713 K4 ,7714 Ks ,7715 ^ 



'50 (3) 

/ ^ / fci,rm fc2,m2 «3,?n3 -k^ ,-m? ^ 



X 



I ki..2 / U3 / U4;5 I ^ ^> 



X (2^^ + 1) 



mi:2 Tn^ J Kms m? —w}) lm4:5 —w? 



€1^2 Jfcl:2 1:2 



X 



1:2 ^ ^3 \ / t4:5 \ „ , ,^m^-fc^ 



Define ^"2 = ip, kP'^ = kp, mP-^ = nip 



/„^,„, 11 -^fca,ma"y - V ^) 11 '-'^<',fe«;W,fea+2W°,m'';^a+2,ma+2 



^ 4+2 ^"+1 



Lemma 29 

P-3 



11 fca,maV ) ii I j„ m"+V 

^ ._-,NS(mi^f-3-fci^P-3) TT 4+2 



a=l a=0 

P-3 
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Proof. Induction 



V- r){«i) r)(^2) n(^3) jj{i4) (_-.^m ( h:2 ^ \( ^ ^3:4 
ki,mi k2,m.2 fcs.ms fc4,m4 v ) \mi-2 —ml \m 7713-4 

mi:4 ^ ■ ^ ^ 



Indeed 

^(.3) ^(..) ^ y (2^^ + l)f^3:4 nff- nz.(f); 



nil + "1-2 + w-s + "1-4 = 



r){^l) r)(^2) r)(^3) r){^4) .-.Nm h:2 ^ \ ( ^ ^3:4 
A:i,mi fc2,m2 ^3,^3 k4,m4 ^ > \rni;2 -mj\m 777-3:4 

/ J / J ki,mi fe2,m2 k 



k'^.rrfi 

4:4 /4:4 / ^1:2 M / ^ 4: 



(2^2 + i)(_i) 

\ ' > 77l3:4 777 y V fe3:4 K J V 77tl:2 — TH, 7 \ 771 7773:4 



,2,^^2 4l^"'3^' -777 ;V"^ "^3:4;^ 

V- n(^i) n(^2) ._-,xm2_fc2 /£3:4 f \ ( ii,2 l\ 

2^^k„rm^k2,m2^k^m^^ ^) U I -mi 

E E <:L.<i°itL (^r, _S) (-1)- 

/>2 ;.2 ™2 V / 



l"^ ,k^ m? ,m\;2 
h:2 n , /^3^4 



^4: 



use m 

Proof, of Lemma [19] We have 



C-.(<0 = En<L.Cum4(C7 

'Tll:p a=l 
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Under isotropy assumption 

p 



Cum„ ( Z\ 



2^ ^ ^) ii I -A;'^ ka+2 k^+^ 

mi;p a=0 ^ ^ 



If isotropy is not assumed, apply 



Gump (Z,7^;-) - (-l)""" ' n (_!„<^ 

p-3 

a=l 

We have 

P-3 



mi:p a=l a=0 
P-3 



a=l 

P-3 



« 4+2 r+1 



a=0 
P-3 

a=l 

see Lemma [29] . ■ 
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